The Cosmology of Ricci- Tensor-Squared Gravity in the Palatini Variational Approach 
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We consider the cosmology of the Ricci-tensor-squared gravity in the Palatini variational approach. 
The gravitational action of standard general relativity is modified by adding a function f{R'^'^Rab) to 
the Einstein-Hilbert action, and the Palatini variation is used to derive the field equations. A general 
method of obtaining the background and first-order covariant and gauge-invariant perturbation 
equations is outlined. As an example, we consider the cosmological constraints on such theories 
arising from the supernova type la and cosmic microwave background observations. We find that 
the best fit to the data is a non-null leading-order correction to Einstein gravity, but the current data 
exhibit no significant preference over the concordance model. The growth of non-relativistic matter 
density perturbations at late times is also analyzed, and we find that a scale-dependent (positive or 
negative) sound-speed-squared term generally appears in the growth equation for small-scale density 
perturbations. We also estimate the observational bound imposed by the matter power spectrum 
for the model with /(R''''R„6) = q(R''''R„6)'' to be roughly |/3| < 0(10"'^) so long as the dark 
matter does not possess compensating anisotropic stresses. 

PACS numbers: 04.50.-|-h, 98.80.Jk, 04.80.Cc 
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I. INTRODUCTION 



The accumulating astronomical evidence for an accel- 
erating cosmic expansion has stimulated many investiga- 
tions into the nature of the dark energy, or other possible 
deviant gravitational effects, which might be responsible 
for this unexpected dynamics (for a review see, e.g., pi]). 
Besides proposing to add some new (and purely theoreti- 
cal) matter species into the energy budget of the universe, 
many investigators have also focused their attentions on 
modifying general relativity (GR) on the largest scales, so 
as to introduce significant modifications in the behaviour 
of gravity at late times when it is comparatively weak. 
One example of the latter sort is provided by the family 
of f{R) gravity models, which had also been considered 
before the discovery of cosmic acceleration (see for exam- 
pie Refs. gii) with reference to alternative forms of 
inflation and the existence of singularities. In Refs. [EQj 
the authors discuss a specific model where the correction 
to GR is a polynomial function of the B?',R°'^Rab and 
R"^'"^'^ Rabcd quadratic curvature invariants (here, R, Rab 
and Rabcd are respectively the Ricci scalar, Ricci tensor 
and Riemann tensor calculated in the standard way from 
the physical metric gab) and showed that there exist late- 
time accelerating attractors in Friedmann cosmological 
solutions to the theory. Barrow and Clifton established 
the general existence conditions for de Sitter, Einstein 
static, Godel universes in theories where the Lagrangian 



is an arbitrary function of these three invariants [3| . 

When the Ricci scalar R in the Einstein-Hilbert action 
is replaced by some general functions of R and R'^'^Rab, 
it becomes necessary to distinguish between two different 
variational approaches to deriving the field equations. In 
the metric approach, as in Refs. [5i,i6|], the metric compo- 
nents gab are the only variational quantities and the field 
equations are generally of fourth-order, which makes the 
theories phenomenologically richer but more stringently 
constrained in most cases. Within the Palatini varia- 
tional approach, on the other hand, we treat the metric 
gab and the connection FJJ^ as independent variables and 
extremize the action with respect to both of them; the 
resulting field equations are second order and easier to 
solve. The Palatini f{R) gravity is also proposed as an 
alternative to dark energy in a series of works [1, [^, . 
There has since been growing interest in these modified 
gravity theories: for the local tests of the Palatini and 
metric f{R) gravity models see [ill [T2 | : f or the cosmolo- 
gies of these two classes of models see |l3l, [H, [H, [H, [13] . 

Both approaches to modifying gravity are far from 
problem-free. In the metric f{R) gravity models, the 
theory is conformally related to standard GR plus a self- 
interacting scalar field [sf , which generally introduces ex- 
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tra forces inconsistent with solar system tests [12|. The 
Palatini approach, on the other hand, generally leads to a 
large (or even negative) sound-speed-squared term in the 
growth equation of the non-relativistic matter perturba- 
tions on small scales. This induces effects on the cosmic 
microwave background (CMB) and the matter cluster- 
ing power-spectra which deviate unacccptably from those 
which are observed [IB, Hi, ITtI . Again, these examples 
highlight the difficulties encountered when trying to make 
modifications to standard GR which are compatible with 
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observations. 



A. The /(R''''R„6) Gravity Model 



In this work we will focus on the Ricci-squared gravity 
models within the Palatini variational approach, which 
we also denote by the /(R^^Rah) gravity. It turns out 
that the Ricci tensor, Rab, and Ricci scalar, R, appear- 
ing in the field equations in the Palatini approach are 
not the ones calculated from the physical metric, gab, 
(we consider the metric gat as the physical one because 
it is this metric which the matter Lagrangian density 
Cm depends on and the energy-momentum conservation 
law holds with respect to) as in GR, and we denote the 
GR equivalents by Rab and R respectively to distinguish 
them from the Palatini quantities). Such a modification 
of gravity has indeed been considered in fl^ and shown 
to give an accelerating cosmology. However, our work 
differs from [3| in that we replace the Ricci scalar in 
the gravitational action with R-|- /(R''''Rab) rather than 
simply /(R''''RQf,), and we concentrate more on the cos- 
mology at the first-order perturbation level, especially 
the late-time cold-dark-matter (CDM) density perturba- 
tion growth. We emphasize the similarity to the Palatini 
/(R) gravity models also. 

Our presentation is organized as follows. In Sec. |TT] 
we briefly introduce the model and outline the meth- 
ods used to derive the background and first-order covari- 
ant and gauge-invariant (CGI) perturbed field equations. 
In Sec. IIIH we present the modified Friedmann equa- 
tion and apply it to a specific family of theories with 
/(R^^'Rah) = a(R"''Ra6)^; the constraints on the pa- 
rameter space (/3, Q„i) from cosmological data are also 
given. Then, in Sec. llVi we analyse the growth of CDM 
density perturbations at late cosmological times for this 
model. Since this analysis shares some similarities with 
the Palatini /(R) gravity models, we also present a simi- 
lar discussion of the latter for comparison. Our discussion 
and conclusions are presented in Sec. |Vl 

Throughout this work our convention is chosen as 
[Va,V6]u'' = R^{; ^lU"^ , Rab = Racb whcrc a, 6, ••• run 
over 0,1,2,3 and c = = 1; the metric signature is 
(-I-, — , — , — ) and the universe is assumed to be spatially 
flat and filled with CDM and black body radiation. 



We will start our discussion with the modified Einstein- 
Hilbert action in the present model. 



S = 



xy/-g 



R- + ./(R-°''R-a6 



2k 



(1) 



in which k — SttGn, with Gn the Newtonian gravita- 
tional constant. Here, Rab = '^abi^bc) given by 



R 



■ab 



^ cd^ ab 



ad cb 



(2) 



and R — g°''''R.ab', note that FJJ^ is a new and independent 
variable with respect to which we extremise the action, 
and is different from the Christoffel symbol F^^ calculated 
using the metric gab- Rab is assumed to be a symmetric 
tensor (if it contains an antisymmetric part then the field 
equation will be spoiled as discussed in [19|) and gab could 
be used to raise or lower its indices. Varying the action 
Eq. ID) with respect to the metric gab (note that ST/Sg = 
as they are independent) gives the modified Einstein 
equations: 



Rn 



1 



2FRSRb,--gab[R+/(R'^''R,b)] = ajT/, (3) 



where F = df/dS with S = K^'^Rab and t/^ is the 
energy-momentum tensor of the fluid matter (CDM and 
radiation). 

On the other hand, varying the action with respect to 
the new variable F^^ with the relation SHab = '^c{SFab)~ 
Db(<5r^^), one arrives at another field equation 



vMetg (5"^ + 2Fg"m,<jg' 



bd\ 



0, 



(4) 



where represents the covariant derivative compati- 
ble to (the covariant derivative compatible to gab 
is denoted, as conventionally, by Vq). Just like in the 
Palatini /(R) models, this equation implies some rela- 
tion between the physical metric gab and the metric gab 
whose Christoffel symbol is F^^. However, because of the 
presence of the second term in the parentheses this rela- 
tion is nontrivial and some further algebra will be needed 
to explicate it. Before doing that, we will present some 
preliminary definitions and expressions, one of which is 
the notation of 3 -|- 1 decomposition. 



II. FIELD EQUATIONS IN /(R^'R^b) GRAVITY 



In this section we briefly introduce the main ingredi- 
ents of /(R°''Rq6) gravity and outline the strategy for 
deriving the general perturbation equations that govern 
the dynamics of small inhomogcncities in the cosmologi- 
cal models that arise in this theory. 



B. The 3 + 1 Decomposition 

The main idea of 3 -M decomposition [13, [HI, [H, 
is to make spacetime splits of physical quantities with re- 
spect to the 4-velocity u° of an observer. The projection 
tensor hab is defined as hab — gab — UaUb and can be used 
to obtain covariant tensors perpendicular to u". For ex- 
ample, the covariant spatial derivative V of a tensor field 
^d - e defined as 

i/'^Tt:^ ^ h-h] ■ ■ ■ hih^, ■ ■ ■ hy^T^:::^. (5) 
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The energy-momentum tensor and covariant derivative 
of the 4-velocity are decomposed respectively as 



Tab = T^ab + '^q{aUb) + pUaUb - phab, 



(6) 



VaUfe = (Tab+^ab+ ^Ohab+UaAb. (7) 

In the above, nab is the projected symmetric trace- 
free (PSTF) anisotropic stress, qa the heat flux vector, 
p the isotropic pressure, aab the PSTF shear tensor, 
•n^ah = V[aM{,], is the vorticity, = V'^Uc = 3a/ a (a is 
defined here as the mean expansion scale factor) the vol- 
ume expansion rate scalar, and Ab = Ub is the fluid ac- 
celeration; the overdots denote time derivatives expressed 
'AS (j> — u°-V a<t), brackets mean antisymmetrisation, and 
parentheses symmetrization. The velocity normalization 
is chosen to be u°'Ua = 1- The quantities TTab^qa, P,P are 
referred to as dynamical quantities and (Tab, '^ab, Q, Aa as 
kinematical quantities. Note that the dynamical quanti- 
ties can be obtained from the energy-momentum tensor 
Tab through the relations 



p 


= TabU'^u'', 




p 






qa 


= h^u'Tcd, 




T^ab 


= hlhfTcd+phab. 


(8) 



Decomposing the Riemann tensor and making use the 
Einstein equations, we could obtain, after linearization, 
the perturbed (constraint and propagation) equations 
[20l. w\. [2^. [23I . Here, we shall not list all of them because 
most are irrelevant for the following discussion; rather we 
will use the linearised Raychaudhuri equation 



^e'-V^Aa + ^{p + 3p) = 0, (9) 



the linearised conservation equations for the energy den- 
sity: 



p+{p+p)e + 1/''qa = 0, 

and the linearised Friedmann equation 
-9^ = np. 



(10) 



(11) 



The above equations are derived and presented for 
standard general relativity, and so the p,p,qa variables 
describe imperfect fluid matter. For general modified 
gravity theories, such as those presented here, the mod- 
ification to GR might be parameterized as an effective 
energy-momentum tensor. In this case the formalism of 
these equations is preserved and one just needs to replace 
p,p,qa by the total effective quantities of the same sort: 



C. The Field Equations in /(R''''R„6) Gravity 

In Eq. (HI), we see that y/det^ {g"'' + 2Fg'"'Ilcd9^'^) is 
a symmetric (2,0) tensor density of weight 1, and so we 
can introduce a new metric gab by means of the following 
relation 



i/det 



gg 



(g'^" + iFg'^'R.dg 



bd\ 



(12) 



where the Levi-Civita connection of the metric gab is just 



rjj^, as we referred to above. 



To go further, we need to express Rob explicitly. This is 
easy to do in principle, because Eq. ([3]) is just an algebraic 
equation for Rab- To see this, let us write the symmetric 
tensor Rab in a general way as 



R„ 



AUaUb + ^hab 



2U(aTb) 



(13) 



where Ua is the 4-velocity of the observer referred to 
above. Substituting Eqs. (l6l[T3l) into Eq. ([3]), we get 

+ 2u(aTb) + Sab 
+ 2F [A'^UaUb + E'^hab + 2(A + S)u(aTb) + 2SEab] 

-i (A + 3S + /) UaUb - i (A + 3S + /) hab 



K [p^Ua 



Ub - p- 



•^ha 



'ab ; ' 



which leads to the following four equations: 

A + 2FA2-i(A + 3S + /) = Kpf, 

S + 2FS2-i(A + 3S + /) = -Kpf 
[1 + 2F(A + S)] Ta = kqI 



(l+4FS)Sab 



'ab' 



(14) 

(15) 

(16) 
(17) 



where f,F are functions of R''''Rab A^ + 3E'\ Thus 
given the specified form of /, and the values of p^,p-^ , 
the quantities A,S can be obtained from Eqs. ([HI [T^ . 
at least numerically. Then, Ta and Sb can also be cal- 
culated from Eqs. (UHl [I7|) provided the values of g£ and 
'^ab given. Note that Ta and Sab are nonzero only at 
first order in perturbation. Taking the time derivatives 
of Eqs. (fT4l [T5| , and using the background values of A, S, 
we could easily obtain A and S by solving the two linear 
algebraic equations^ Similarly, VaA and VaS could be 
worked out (here, V is the spatial derivative). In what 
follows, we shall assume that A, S and their derivatives 
have been calculated. 

The next step is to find out the relation between gab 
and gab. We could rewrite Eq. as 



Vd^g'^" = y/d^g^'i6'^ + 2FRl). 



(18) 



Taking the determinants of both sides and equating we 
get 



detg = detg-detP, 



(19) 
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with 



5b 



2FRg. 



(20) 



Thus, we conclude from Eq. (|18|) that 



ab ^ V3et^ acpb 
^ /J — I — -J c 



Vdetj 
1 

VdetP 



Vdetg 



(21) 



(22) 



Obviously, detP and P ^ need to be evaluated respec- 
tively. For detP, we have 



detP 

det[(l + 2FA)u''ub 



(1 + 2FS)/ig 



To calculate this, let us write goo — a^(l + 2'^),goa — 
gaQ = a?Ba,gai3 = [(1 + 2HLha0 + '^Hxap] and uq = 
a(l + ^'), Met = —a{va — Ba) where a, (3 run over 1, 2, 3, 
, Ba, Hl and HxafS are first order metric variables of 
which Hxap is traceless, is the spatial component of 
Ua, and "fafS is the metric of 3 dimensional flat space. As 



a result 



,00 



(1 

fa/3 



2^0,^ gap ^ 



{l~2HL)-i°'f^~2H^^ and 1*° = a-i(l-*),"" = -a"^w". 
From these expressions the components of /ijj can also be 
obtained and one can substitute all these quantities into 
the above equation to get detP. Since and are 
only of first order and because Safe is traceless, it is then 
not difhcult to see that up to first order (note that the 
facts M°Ta = MaT° = and u°Sab = indicate that 
TO = To = and 2° = Eg = 0) 



detP = (1 + 2FA)(1 + 2FS 



(23) 



For (P~^)^, we know that it is symmetric as the inverse 
matrix of a symmetric matrix, and so could be written 
as 



(P-1)° = Au''ub + Bhl + 2u^''C^+Dl. (24) 



Using 



= (1 + 2FA)u''mc + (1 + 2F'~)hl + AFu^^T^) + 2FY.\ 



and 



it is then easy to obtain, to first order, that 
A 
B 

Ca 
Dab 



l + 2i^A' 
1 

1 + 2FS' 



2F 



(l + 2i^A)(l + 2FS) 
2F 



T • 



(1 + 2F^f 



„x9 E(jb 



(25) 



As a result, we have now the relations between the two 
metrics gab and gab and their inverses as 



gafc — ^9ab + £,ab 



Qab 



where 
A 

iab 
^ab 



v/(l + 2FA)(l + 2FS), 
1 + 2FS 
1 + 2FA' 

X{UJ - l)UaUb - 4VwFU(aTh) 



2F . 



1/1 \ r, b 1 2F 



(26) 
(27) 

(28) 
(29) 

(30) 
(31) 



A discussion of how the two Ricci tensors Rab(rg^) 
and Rabidab) are related to one another is given in the 
appendix, with the help of which the Einstein equation 
Eq. ^ can be rewritten as 



Rab 



:gabR 



ab 



kt: 



where 



^9abif + SR) 



SRab 



ab 



2FR^R, 



cbi 



(32) 



(33) 



and (see the appendix for a definition of the tensor 7^^) 

SRab = yci:b-^bi:,+ltblcci-ltclbd, (34) 



SR = g^'^SRab 



(35) 



With the aid of Eqs. (0 [21 [Ml [S3 one could iden- 
tify p^f-^,p^f-^,q^f-^ and tt^I^ and express them in terms 
of oj, X, F, A and S, which are functions of p-^ and 
(c.f., Eqs. IHl [13 [Ml [29])), and Ta,Eab which are also 
functions of qI^t^Ii, (c.f., Eqs. ^JM \T7\ ). However, from 
Eqs. ([26] - [311 [Ml EOl) one can see that this process will 
involve a lot of calculation. In the present work we will 
not perform detailed numerical calculations of the per- 
turbation equations of the /(R''''Rq6) model; instead, in 
the next two sections of the paper we will: 

1. Study the background evolutions of general Ricci- 
squared gravity models. As an example, we will con- 
sider a specific family of theories with /(R''''Raf,) = 



5 



a(R°''Rab)'^, and constrain the allowed (a,/3) parame- 
ter space using data sets supernovae (SNe) luminosity 
distances and the CMB shift parameter. 

2. Present a simple argument to show that this class of 
modified gravity theory, like those arising in the Palatini 
/(R) theory, generally possesses a scale-dependent effec- 
tive sound-speed-squared term which affects the growth 
of CDM density perturbations and thus influences the 
matter power spectrum [l^, [l6|, |T3| on small scales. 



III. THE COSMOLOGICAL BACKGROUND 
EVOLUTION 

In order to analyse the cosmological background evolu- 
tion we can neglect the and tt^j, terms, and hence the 
quantities Ta, Sab- As a result, the equations are greatly 
simplified. 

We are interested in the modified Friedmann equation 
in the present model. From Eqs. ([8l fTl ] [33 |) we have 



up 



(36) 



in which H = ^9 is the Hubble expansion rate and p 
is expressed as 



up 



Kp 



f 



i/ - 2FA'' + ^SR - 6Rabu''u\ (37) 



with SR and 6Rab given in Eqs. p4l[35|) . After a lengthy 
calculation, and using Eq. Q to eliminate the term 9 + 
3A/2A which appears in Eq. fST]) . we obtain the following 
simple result 



H 2 



g(A-3^S). 



(38) 



There are two interesting points regarding Eq. (|38p . 
Firstly, we see that only w, and not its time derivatives 
uj or a), enter the equation. Secondly, the second-order 
derivative of A does not appear either; to see the con- 
sequence of this, note that since A = d\{p^) / dp^ p^ = 
—<^d\{p^)/dp^p^H (with = 3 for matter and <; — A 
for radiation), the dependence of A on pf could be ex- 
pressed in terms of p-^, e.g., in radiation-dominated era 
pf = P'^/3 and in matter-dominated era p^ — 0, and so 
we have A cx 77. Consequently, Eq. ((38)) has the form 



e(p^/), 



(39) 



where is a complicated function of p^ and p-^ (at late 
times p^ ^ and it becomes a function of p^ alone) . 

As we discussed above, knowing pf and p-^ means that 
we know A, S, A and lu. Thus, given a specific form 
for /(R°^Rat), Eq. ([55)1 completely determines the back- 
ground cosmological evolution of the model. As a partic- 
ular example let us consider the case of 




TOOO 



FIG. 1: (color online) The time evolution of the effective 
equation of state for the /(R''''R„b) = a(R"''Rab)'' model. 
For clearness we only show the cases of P — 0.2 (the dashed 
curve) and (ACDM, the solid curve). The other parameters 
defining the cosmology are Qrn = 0.27 and ilr = 8.5 x 10~^. 
Clearly an accelerating phase following the standard matter- 
dominated era can be realized within this model. Note that 
the redshift range 1 + z < 1 corresponds to the future of the 
present epoch, where the Universe evolves into a de Sitter 
stage. 



where a and f3 are the model parameters. Note that here 
R^^Rafc = -I- 32^ is always non-negative, and /3 = 
corresponds to a picking the standard ACDM cosmology 
of GR. 

For convenience, we shall define the following dimen- 
sionless quantities 



/ 
A 



F 



J_ 
A 



npr 



4/3-2 



/(R'^'Rab) = aiR'^'Rab) 



(40) 



then Eqs. [T5|) could be rewritten as 
A + 2FA^ - i (a + 3S + = 3fi,; 
S + 2FS2 - i (a -f 3S -)- /I = -fir 



(41) 

-3a-, (42) 
(43) 
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where 



62 



and than Eq. 



1 



reduces to 



2A 



IP 

HI 



(44) 
(45) 

(46) 



where Ap/ = dX/dp^ . 

In this paper we will set fir = 8.5 x 10^^, so today we 
have H^/H^ = 1 and there are 3 equations (Eqs. 

for the 5 parameters A,S,a,/3 and f7,„. Therefore, 
we are able to express all the other quantities in terms 
of f3 and flm, which can therefore be treated as the two 
independent degrees of freedom of our model. Note that 
a is a constant, and once evaluated at the present day, it 
could be used all through the cosmic history, which helps 
determine A,S at arbitrary times. 

In Figure [1] we have plotted the effective equation of 
state, defined by Weff = — 1 — |-^ = — 1 — (where 
a star-superscript denotes the derivative with respect to 
log(a)), as a function of the redshift. The values of f3 
are indicated beside the curves. At early times, when 
the f {R°'^Rab) corrections are negligible, the models all 
mimic the ACDM evolution of Wef / , and the same thing 
happens in the future. This is because during this era the 
matter (relativistic and non-relativistic) is greatly diluted 
so that the right-hand sides of Eqs. ([1^ both vanish; 
consequently, we can solve them to show that A = S = 
const, and so f{R'^^Rab) is also constant. The deviation 
from ACDM occurs mainly at intermediate times, that 
is, in the recent past and future. 

We now use the observational data on the background 
cosmology to constrain the parameter space (in the /3 — 
rim plane) of the present model. For this we jointly use 
the 157 measurements on SNe luminosity distance in the 
Gold data sets of Riess et al. [H] and the CMB shift (R) 
parameter. The SNe luminosity distance is expressed as 



dL{z) = (1 + 2) 



1 + z 



Ho 



du 
du 



(47) 



where E{z) = H{z)/Ho. The measurements supply the 
extinction-corrected distance modulus fiQ = 51gc?L -I- 25 
(with in units of Mega-parsecs) and its uncertainty, a, 
for individual SNe, so that the standard minimization, 
defined by 

X - 2^ -2 (.48) 

is easy to implement, where fj,p is the theoretically pre- 
dicted distance modulus. As Hq appears only as it does 



0.40r 




FIG. 2: The constraints on the parameter space of flm and /3 
in the present model from joint SNe and CMB shift parameter 
data sets. The grey and light grey regions represent the 68.3% 
and 95.4% confidence contours respectively. The white circle 
{Q,m = 0.265, /3 = 0.07) is the best-fitting parameter of our 
model, and the star is the concordance ACDM model. 



in Eq. (|47p . we could marginalize over it by integrating 
the probability density p(x) cx; exp(— x^/2) for all values 
of Hq. For the CMB R-parameter, defined as 



R 



dz 



(49) 



we adopt the observational value R°^^ = 1.70 ± 0.03 at 
Zdec ~ 1089 from [2^. Note that this does not depend 
on the specified value of Hq. 

Our constraining result is shown in Figure [21 where 
we have shown the 68% and 95% confidence regions re- 
spectively. The constrained intervals are roughly 0.20 < 
n,n < 0.36 and -0.13 < (3 < 0.29 at the 95% confi- 
dence level, with the best fitting values being (/3, flm) — 
(0.07,0.265) with x^/dof ~ 1.126. Also note that the 
concordance ACDM model (the white star) lies within 
the 68% confidence region of our constraints. 

Thus, we see that the background cosmological data 
is able to constrain \(3\ to be of order 0.1. In the next 
section we will bricfiy investigate the possible constraint 
from the growth of dark-matter density perturbations, 
and show that this may provide a potentially more strin- 
gent restriction on (3. However, considering that this lat- 
ter limit depends on the properties of the dark matter, 
our background constraints given in this section are less 
model dependent. 
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IV. EFFECTS ON LATE-TIME CDM DENSITY 
PERTURBATION GROWTH 

In this section we study the effects of the /(R^^Rafc) 
corrections to GR on the CDM density perturbation 
growth. We start by recalling the case of /(R) gravity 
because it shares some similarities with the /(R^^Rat) 
one, while being technically simpler than the latter, and 
because a similar analysis for the former is still missing 
from Refs. (see however [15] for a slightly different 

treatment). 



can be rewritten as 

Rab 



in which the effective total energy-momentum tensor is 
given by 



rritot 
ab 



-kTI 



A 

Using Eq. 



^5„fc(R + /) - i.gab(R - 6R) - 5Rab.{53) 
, we can now identify 



A. The Case of Palatini /(R) Gravity 

Recall that in our simplified model the universe is filled 
with CDM and radiation, and at later times the radiation 



1 



/ - 



PCDM- 



energy density is negligible, so p 

Taking the spatial derivative of the Raychaudhuri 
equation Eq. ([9]) (with the p,p there being replaced by 



p'^-'-' ), and working in the 
CDM frame (where the observer is comoving with CDM 
particles and thus A = 0) [2^, we have 



A" 

^CDM 



HA'r 



CDM 



0, (50) 



where Acdm is the CDM density perturbation contrast 
that is defined through VaPCDM = PCDuJ^k^^Qa^ 
and Ti = |a is the Hubble expansion rate with respect 
to conformal time (note that a prime denotes the con- 
formal time derivative, and a dot the cosmic comov- 
ing proper-time derivative); X^X^ are respectively the 
harmonic expansion coefhcients for VaP and V^p (de- 
fined via V.p = Efe l^Ql and V.p = l^^Ql 
[26l|). Clearly we need to know about X^^i and XT^ ^^^ 
which arise from the /(R"''Rafc) modifications to GR 
(c.f. Eq. (33)). 

In the Palatini /(R) model, in which the Ricci scalar 
R in the gravitational action is replaced with R -|- /(R) , 
Eqs. (El HZ]) stiU hold, but with (see for example ie|) 



A = 1 + 7^ 



9R' 



Cab 
Cab 



0, 
0. 



(51) 



Then, with the help of the calculations in the appendix, 
it is straightforward to show that 

Rob = Rab + ^Rab 

3 11 

= -Rab + ^VaAVbA--VaVbA-—5abnA (52) 

where □ = V^, and the modified Einstein equation (l6j . 



ARab- 25ab(R + /) = ^Tl^. 



Kp 



f 



-Kp 

1 

2 
3A2 
2A2 

/ 



1 

2A 



R+^QA-^V^AVaA 



-Kp 
1 

'2 



A 

A ' 
1 

2a' 



R + /) 



(54) 



Thus 



k(p*°* + 3p*°*) = ^A + H^A 
A A 



(55) 



(56) 



in which • • • represent the terms not involving second 
order (time and spatial) derivatives. 

The reason why we keep only two second derivative 
terms explicitly on the right-hand side of Eq. ([56]) is that, 
after taking the spatial covariant derivative, the first term 
contributes a A^^jy; piece to Eq. ([50]) while the second 
term contributes a fc^AcDM piece. None of the remaining 
terms in • • • contribute these two pieces to Eq. ((50|) . To 
be more explicit, recall that A = A(/3cdm) in the model, 
so 



VaA 



dXjpcBM) 
dpcDM 

X 



VaPCDM 



PCDM 



-Va/CCDM 



A 



'SHpcBU 



VaPCDM 



A \ ^ k ^ 



(57) 



where we have used Eq. (fTU]) to background order. As a 
result, 



k{X*°^ 
X 



XH 



A" 



X 



CDM 



3AiJ 



k Acdm 
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and Eq. ((50l) can be recast into the form 
A 



1 



2\H 



^CDM + [• • • ]^CDM 



A 



6AiJ 



A 



CDM 



which, after rearrangement, gives 



A" 

^CDM 

A 



3(2AiJ + A) 



CDM 



iCDM 



(58) 



Here, • • • denotes comphcated terms that are determined 
completely by the background evolutions of the model, 
and are unimportant for our analysis here. What is es- 
sential in Eq. ([55]) is that it tells us that, as long as the 
quantity A does not vanish, in general there will appear 
an effective sound-speed-squared term for the growth of 
matter density perturbations. Depending on the sign of 
A, this sound-speed-squared term could be either positive 
or negative, in both cases the small-scale density pertur- 
bation growth becomes extremely scale-dependent, alter- 
ing the shape of the matter power spectra significantly 
[Tsl . [TgI . ITtI . Notice that the terms in [• • • ] could also 
modify the evolution of density contrasts differently as 
compared with the prediction in standard general rela- 
tivity, but in a scale-independent manner, and at small 
scales their effects are subdominant. 

One more comment on the modified Friedmann equa- 
tion in the Palatini gravity models is appropriate. Using 
Eq. (fTTj) with p replaced by the given in Eq. ([54]) , we 
can see that only the A0, A'^ terms are involved and the A 
terms cancel (we only consider terms to background order 
here). Since A = d\{p^)/dp^ — ~(^d\{p^) / dp^ p^ H we 
see that A0, A^ cx and could be moved to the left-hand 
side of Eq. (11); the remaining terms on the right-hand 
side are also functions of p^ only and so Eq. ([55]) is also 
realized. To be more explicit, the modified Friedmann 
equation in Palatini /(R) gravity is [l^ 



H 



2A 



1 

6A 



[^^ipf +3pf)-{-R + f)] 



which can be shown to be just Eq. ([55]) ii lu — 1 there, as 
expected, because the metrics then take the same form 
(of course, the definitions of A are different in the two 
cases). 



simpler theory described in the last subsection shows that 
what is most relevant for our analysis are the second or- 
der (time and space) derivative terms (note that the term 
9, if exists, can also contribute to Acdm because Va^ 
contains Z where Z is the Harmonic expansion coefficient 
of Va^ via VaO — ^ZQI, and because aAcDM = —kZ; 
however it turns out no 9 terms appears in K(p*°*-|-3p*°*)), 
which are straightforward to identify. 

We shall formally repeat the procedure of the last sub- 
section. Note that the quantities Ta, Sab are determined 
by ql^ , TT^j, which in our case are due to the radiation mat- 
ter species. At late times the radiation energy density is 
negligible so that, to a good approximation, qIjTtIij and 
thus Ta,Sa6, vanish. As a result, Eqs. ((26l [27|) become 



A.gab + X{UJ - l)UaUb, 



1 /I 

aU-^ 



(59) 
(60) 



with A, LU defined in Eqs. (^511^5)1 . Meanwhile, since p-^ = 
PcDMi we have p-'' = at late times; F, A, S , and hence A 
and u), become functions of pcDM only, i.e., A = A(pcdm), 
LU — w(pcdm)- The analysis in Eq. (|57p then also applies 
to A and oj here, so 



k 

v.(Ac.) = -^^^^AcDMQ^ 

6H — ' a 



Now, from Eq. ((33l) we obtain 



«;(p*°* + 3p*°*) 

2k/9*°* - k(p*°* - 3p*°*) 

k(/ -f 3/) - 2F{A'^ - 



3S2) - / - 2SRabu''u\61) 



where the relation RJjRch = A'^UaUb + ^^hab and Eq. ^ 
are used; 6Rab is given in Eq. ([M]) . Note that SR does 
not appear in this formula, and what we need to evaluate 
is just the collection of second-order derivative terms in 
5RabU°'u^ ■ After some manipulation we obtain a similar 
expression to Eq. ([SS)) for the Palatini /(R) model: 



.3p'°*) = |A+iv2(Ac^) + --- (62) 
A A 



The following analysis then completely parallels that for 
the Palatini /(R) model and the CDM density perturba- 
tion growth equation can be shown to be (like Eq. ([58]) ) 



B. The Case of Palatini /(R"''R„(,) Gravity 



Now consider the /(R" Rq^) gravity model. As dis- 
cussed in the last section, the detailed forms of p*^^^ .p*^^^ 
are very complicated, but fortunately we need not eval- 
uate the full formulae explicitly. Our experience of the 



A'cDM + [■ ■ ■ ] ^CDM 



[Xlu)- 



3{2XH + X) 



iCDM 



(63) 



Again the [ • • • ] in Eq. denotes the terms which are 
completely determined by the background evolution and 
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are not of interests to us here. Thus we see that, simi- 
lar to the case of Palatini /(R) gravity model, a scale- 
dependent sound-speed-squared term also appears in the 
Palatini /(R''''Rah) gravity model, whose sign depends 
on (Aw)'. Note that in the case of w = 1 the metric gab 
(Eq. (1591 160]) ) has the same form as that in the Pala- 
tini /(R) model, and Eq. (63) reduces to Eq. (58) as 
expected. 

In the general /(R^^Rat) gravity models it is possi- 
ble that (Aw)' 7^ 0, thus the scale-dependent effective 
sound-speed-squared term influences the matter pertur- 
bation growth and could alter the shape of the predicted 
matter power spectrum. Our previous knowledge derived 
from refs. [IB, O suggests that this effect might al- 
low observational data to place stringent constraints on 
the parameter space of these theories. In fact, we can 
give a rough estimate of how stringent the constraint 
can be. Consider first Eq. ((55)) for the Palatini /(R) 
model: since on small scales the fc^ term dominates the 
other terms in front of Acdm (the ' • • •' terms), it is ob- 
vious that the magnitude of the quantity A/3(2AiJ -I- A) 
determines the deviation from ACDM results. The ob- 
servational constraint on the model parameter /3 (recall 
that /(R) = a(-R)'^), as shown in Refs. [iB, [3, is 
< 0(10-6 ^ 10-5); we take ~ IQ-^ and = 0.3 
for illustrative purposes, and find that |A/3(2A7?-|- A)| ~ 
0(10-"^ - 10-6) for the relevant redshift range of < 
z < 10. In the analogous case of Palatini /(R^^Rafc) 
gravity, Eq. (p5)) . we have (Ati;) /3(2AiJ-l- A) as the domi- 
nant term instead of X/3{2XH + X), and so the constraint 
should be |(Aw) /3(2AiJ + A)| < 0(10"^ - IQ-^) in the 
same redshift range. Again, taking flm = 0.3, for the 
theories with /(R'^^'R^b) = a(R"''Rah)'^ we find that to 
satisfy the above constraint must also be limited to 
be 0(10-5) or even smaller. 

Thus, we conclude that the Palatini /(R^^Rat) models 
may be constrained by the observational data on the mat- 
ter power spectrum just as stringently as are the Pala- 
tini /(R) models. Yet, we should note that a more exact 
quantitative constraint can only be obtained by exploit- 
ing a full parameter-space search as done in T6], and that 
the above conclusion depends on the assumption that the 
CDM particles have vanishing anisotropic stress. If, in 
contrast, the dark matter particles admit an anisotropic 
stress, in a manner similar to that prescribed in [27[, 
then the effective sound-speed-squared terms might be 
canceled and leave no significant traces. 



V. DISCUSSION AND CONCLUSIONS 

We considered a general class of modified gravity mod- 
els where the Ricci scalar in the gravitational action of 
GR is replaced by a function R -I- /(R^^Rafc) and the 
field equations are derived using the Palatini variational 
approach, i.e., treating the metric gab and connection 
as independent variables so that the action is varied 



with respect to both of them. The strategy for deriving 
the cosmological equations at both the background (zero- 
order) and the first-order perturbation levels is outlined. 
The main step in this process is to determine the metric, 
gab, whose Levi-Civita connection is F^^, relate it to the 
physical metric gab and thereby fix the relation between 
Hab and Rab- Then, the correction to GR is treated as 
a new effective energy-momentum tensor while the field 
equations take the same form as in GR. The formulae 
laid down here might be useful for the numerical imple- 
mentations of such modified gravity models. 

We also investigated in detail a power-law correction to 
the usual Einstein-Hilbert action given by /(R^^'Raf,) = 
a(R"''Rab)'^. We used the SNIa luminosity distance and 
CMB shift parameter data to constrain its (independent) 
two-parameter space (/3,J7m)- We found that at 95% 
confidence level /3 € [-0.13,0.29] and fl^ e [0.20,0.36]. 
A slightly positive value of /3 (/3 ~ 0.07) is preferred 
by the data used. However, the standard ACDM model 
(equivalent to (3 — 0) with fl„i = 0.27 is still within 
the 68% confidence contour. Hence, although the best 
fit to the data is a non-null leading-order correction to 
Einstein gravity, the current data exhibits no significant 
preference over the concordance ACDM model of GR. 

The late-time growth of matter density perturbations 
in general Palatini /(R'^''Rafc) gravity models was also 
studied. It was shown that the equations governing this 
class of models look very similar to that in the Pala- 
tini /(R) models. In particular, there exists a scale- 
dependent effective sound-speed-squared term in the per- 
turbation growth equation which may be either positive 
or negative, depending on the background evolution in 
both models. In the /(R) case it is well known that these 
terms can lead to strong scale-dependence of the matter 
power spectrum, which is highly constrained by observa- 
tional data p^ . [l6| , and we expect a similar feature to 
exist in the /(R'^^'Rob) case. We estimate that this will 
produce a strong observational bound of |/3| < 0(10^^) 
unless some exotic properties are added to the dark mat- 
ter candidate [27| . 

As the final remark, we give a brief comment about the 
static and spherically-symmetric solutions of the present 
model. The analogue for the Palatini /(R) model was 
considered in [30*1 and the authors found that the exte- 
rior spherically-symmetric vacuum solutions are unique. 
Here we just want to point out the /(R°''Rab) model 
also shares this feature. In fact, in the vacuum where 
pf — pf = it is easy to show that for our model 
A = S = const, are uniquely determined by Eqs. (|141 
[TS)) and so is Rab. The full consideration of a static sys- 
tem also needs the interior solution and its matching the 
exterior, which is beyond the scope of this paper and will 
be left for further investigation. 



Appendix 

In this Appendix we present the relation between 
H-a6(r^^) and Rab{9ab) if the two metrics g^b and gab 
satisfy the following relations 

gab = A^ab + ^afc, (64) 
™(ib ^ ab I /-ab 

g = 1^5 +C (65) 

where A is a scalar function and ^afc , Cafe symmetric ten- 
sors. 

Firstly, the requirement 

§•^^806 = g'^'gcb = <5? (66) 

implies that 

Kab + \u + CaU = 0. (67) 

Then, with some algebra, and using Eq. (|67p . we can 
easily show that (here a comma denotes the ordinary 
derivative) 

nc(gafc) 
— 2^'^'^ {&bd,c + ^cd,b + ^bc,d) 

[(AffM + 6<i),c + (Aged + ^cd),b + (A.g6c + 6c), d] 
= rgc(5afc)+76c, (68) 
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where the difference between F^^ and Ffec", denoted 7^^, 
is defined by 

76c ^ ^KV.A + <5,"VbA-5bcV''A] 

+ ^ [Ct'VcA + CVbA - 5bcC'*VdA] 

+ ^C'"'[Vc6d + V66d- Vd6c] (69) 

and is a true tensor, as expected. From the definition in 
Eq. ^ we can thus derive that 

Rafc - Rab + ^clU - ^bllc + <blcd - 7ac7M (70) 

Note that Rab differs from Rab by a rank-2 symmetric 
tensor. 
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